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Abstract Zero modes of rotationally symmetric vortices in a hierarchy of 
generalized Abelian Higgs models are studied. Under the finite-energy and 
the smoothness condition, it is shown, that in all models, n self-dual vortices 
superimposed at the origin have In modes. The relevance of these modes 
for vortex scattering is discussed, first in the context of the slow-motion 
approximation. Then a corresponding Cauchy problem for an all head-on 
collision of n vortices is formulated. It is shown that the solution of this 
Cauchy problem has a - symmetry. 
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1. Introduction 

Since their discovery vortices in the Abelian Higgs model have at- 
tracted much attention. This is mainly due to the fact that the static so- 
lutions of this model describe flux tubes in superconductors. As objects in 
2-dimensional space, they also provide simple examples of topologically non- 
trivial structures. Recently , the Abelian Higgs model was generalized, and 
rotationally symmetric vortices were found in these models [[J. The gener- 
alised Abelian Higgs models form a hierarchy of 2 dimensional U(l) models 
lablled by an integer parameter p. The significance of p is that the p-th 
member of this hierarchy has been derived by subjecting the p-th member of 
the hierarchy of 4p dimensional scale-invariant SO(4p) Yang-Mills models to 
dimensional descent 0. The p — 1 members of both these hierarchies are the 
usual Abelian Higgs model and the Yang-mills model in 2 and 4 dimensions 
respectively. 

A mathematically rigorous proof for the existence and uniqueness of 
the rotationally symmetric self-dual vortices of the hierarchy of generalised 
Abelian Higgs models was subsequently given in || . These rotationally sym- 
metric solutions describe vortices superimposed at the origin, and provide 
us with objects qualitatively similar to, but quantitatively different from 
the vortices in the Abelian Higgs model. In the case of the Abelian Higgs 
model, actually, a 2n-parameter family of vortices exists 0. In this paper, 
we show that, at least near the rotationally symmetric vortices, there is also 
a 2n-parameter family of generalized vortices. 

In the Abelian Higgs model, the knowledge gained from the study of the 
zero modes has been used to investigate the scattering of vortices. Vortex 
scattering has been studied in the context of the slow-motion approximation 
[I], and with the help of the theory of partial differential equations @. Among 
the most interesting processes studied is 90° scattering of 2 vortices in a head- 
on collision (or its generalization: - scattering in all head-on collisions of n 
vortices.) We will show that some of the results, in particular the - symmetry 
of the scattering process, generalize to generalized vortices. 

The paper is organized as follows. In Sec. 2, we highlight those aspects 
of the generalized Abelian Higgs models, and of their rotationally symmetric 
vortices, important to our discussion of the zero modes. In Sec. 3, following 
steps taken by Weinberg in the case of the Abelian Higgs model |8j , we derive 
the fluctuation equations for rotationally symmetric generalized vortices. We 
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show that these equations have a 2n-parameter family of smooth finite-energy 
solutions. In Sec. 4, the scattering of n generalized vortices shortly before and 
after they form a rotationally symmetric vortex is discussed, first, briefly, in 
the context of the slow-motion approximation. Then a corresponding Cauchy 
problem is formulated. It is shown that its solution has a - symmetry. 

2. The Models and their Radially Symmetric Solutions 

The models we study are the generalized Abelian Higgs models in (2+1)- 
dimensional space-time, given by the Lagrangian densities, 



for any integer p > 1. For p = 1, (1) reduces to the usual Abelian Higgs 
model. is the complex Higgs field, D^cj) = <9 M + lA^cfi, /i = 0,1,2, is 
the covariant derivative, and the gauge fields F^ u are defined in terms of 
the real gauge potentials as F^ u = d^A u — d u A^ = d[^A u ^ /i, v = 0,1,2. 
The square brackets mean antisymmetrization in the indices. For any tensor 
Tftv..., (T^ u ) 2 means T^ T^---. The indices are raised and lowered with the 
metric tensor g = diag(—l, +1, +1). Our task is to find, and study, smooth 
finite-energy solutions of the corresponding Euler-Lagrange equations. 

As in the case of the Abelian Higgs model, the Euler-Lagrange equations 
of the generalized Abelian Higgs models are solved by certain first-order 
equations. For the Abelian Higgs model, these equations were found by 
Bogomol'nyi ||. For the Lagrangian (1), the first-order equations arise as 
follows: In the case A = 0, and for time-independent Higgs field and gauge 
potentials, the Lagragian can be written in the form, C^) = j 2 _|_ J 2 + e ij ^f2 J , 
where J 2 and J 2 are the positive definite terms, 



£(P) 



{rj 2 _ |0|2)2 ( P -2) ([(?? 2 _ |0|2)F ^ + <p _ 1)D[ ^ Dj/](j)] 2 

+Ap(2p - l){rf - |0| 2 ) 2 |ZV| 2 + 2(2p - 1) V - |0| 2 ) 4 ), 



(1) 



J 2 



(^ 2 -|0| 2 ) 2(p - 2) (((^-|0| 2 )^ 
+i{p - l)D^D j]( l>) - (2p - l)e t ,( V 2 - |0| 2 

2p(2p - 1)( V 2 - \<P\ 2 f^\D l( p - ie l3 D^\\ 



2\2\2 



(2) 



(3) 



and where, 
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Vl j = A(2p-l)ri^ 2p -^A j 

(4 ) 

The special structure of the Lagrangian, is no accident, but arises 
naturally through dimensional reduction of the generalized Yang-Mills model 
on R 2 x S ip ~ 2 0. The divergence, e^c^f^, is the dimensionally reduced form 
of the Chern-Pontryagin density. 1 = and J = stem from the self-duality 
equations of the generalized Yang-Mills theory [§]. Setting I and J equal to 
zero, means minimizing the energy in a given tological sector, and yields the 
desired first-order equations. We have simplified the task of finding time- 
independent solutions with A = 0, to solving the following equations: 

( v 2 -\^\ 2 )F tJ +t(p- \)D^D^<j> = (2p-l) e ^ 2 -|0| 2 ) 2 , (5) 

A0 = leijD^. (6) 

For the rotationally symmetric ansatz, 

<fi = r)g(r) exp(— (7) 
jd{r) — n 

A% = tijX — , (8) 

with integer vortex number n, Eqs (5) and (6) reduce to, 

(1-/)^ = ^(p-l)aV-^-l)r(l- 9 2 ) 2 , (9) 

d 4- = 2£. (io) 

dr r 

In Ref. 3, for all n, existence and uniqueness of a solution to these equa- 
tions with the desired asymptotic behaviour was shown. For small r, a goes 
to n, and g goes like C n r n , where the C n are constants which have to be 
determined numerically. For large r, g goes to 1, and a/(l — g 2 ) goes like 



(2p — 1)/ (2p)rjr. In the next section, we will study the zero modes of these 
rotationally symmetric solutions. 
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3. The Fluctuation Equations and the Zero Modes 

We write the fluctuations about the rotationally symmetric solutions in 
the following form: 

5<p = rjg(r)h(r,9) exp(—m9), (11) 

5A 1 = -(-b(r, 9) sin 9 + c(r,9) cos 9), (12) 
r 

cos 9 + c(r, 9) sin 9). (13) 

Here, b and c are real functions, and h = h 1 + ?/i 2 is a complex function. To 
count the number of modes, we will later fix the gauge by setting h 2 equal to 
zero. To discuss the smoothness of the modes, we will have to use the gauge 
freedom and find suitable functions h 2 . 

Equations (5) and (6), linearized in the functions h, b and c, then yield 
the fluctuation equations. Equation (6) leads to expressions for the functions 
b and c in terms of the function h: 

dh 1 dh 2 

Or 89' K V 

dh 1 dh 2 

Using these expressions, we obtain from eq. (5) the following equation for h: 
1 d . dh\ 1 d 2 h} A(p-l)ag 2 dh} 2/n , x2 4(p-l)a 2 Nrl 

(16) 

Note that h 2 does not occur in Eq. (16). Therefore, the finite-energy solutions 
of Eq. (16), which lead to smooth fields in a suitable gauge, yield the zero 
modes. There are no other modes, since we can impose the gauge condition, 
h 2 = 0. 

We attempt to solve Eq. (16) in terms of a Fourier series for h 1 , 

oo 

h\r,9) = Y / (h£\r) cos k9 + h { *\r) sin k9). (17) 

k=0 



5 



The Fourier coefficient functions, have to satisfy the equation, 

(Phf .1 4(p - l)ag 2 .dhf ,k 2 n(n 1N 2 2 4(p-l)oV uW n 
ar 2 r r(l — g z ) ar r z r^(l — (?) 

(18) 

We want to find all solutions such that, in a suitable gauge, the rotationally 
symmetric background fields plus the fluctuations are C°° functions on R 2 
with finite energy. 

For small r, Eq. (18) reduces to, r 2 y" + ry' = k 2 y, where y = h®. 
The general solution to this equation is, y = C\r k + C2r~ k for k > 0, and 
y = c\ + C2 lnr for k = 0. The uniqueness of radially symmetric vortices rules 
out nontrivial fluctuations for k — 0. (We assume here that the uniqueness 
result, rigorously proven in Ref. 4 for p = 2, is in fact true for all p.) For 
large r, 

V p p 

holds, with general solution, 



y = c 3 exp(W-(2p- 1)1 rpr) + c 4 exp(-W^^ — —ryr)- ( 20 ) 

V v V v 

Because of the finite-energy condition, we have to set C3 = 0. C4 is an 
arbitrary coefficient. One can now show that all exponentially decreasing 
solutions at infinity lead to solutions at the origin with C2 7^ 0. Assume 
the contrary: Then for C4 > there must be a maximum with positive y. 
However, Eq. (18) shows that, when the first derivative vanishes, the second 
derivative is positive. This is a contradiction. (An analogous argument holds 
when C4 is negative and y has to attain a minimum.) Hence, acceptable 
nontrivial behaviour at infinity leads to an r~ k term at the origin. For k > n 
this implies that the energy is infinite. 

We are left with the In modes (with all the other fif \ for j ^ % 
or I ^ k, equal to zero) for < k < n and % — 1, 2. Supplemented with 
suitable functions h 2 , these modes lead to the following smooth finite-energy 
solutions: For the first set of n functions, the fluctuations are of the form, 

5<p = rig(r)h£\r) exp(-i(n - k)6), (21) 
5A 1 = -(^- + ^ ) )sin((k-l)6), (22) 
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5A 2 = -(— k - + cos((£; - 1)0). (23) 

ar r 

For the second set of n functions, the fluctuations are of the form, 

5(f) = —iT)g(r)hjp (r) exp(—t(n — k)9), (24) 

J!, ( 2 ) t, 

M 1 = (^ + -4 2) )cos((A;-l)0), (25) 
ar r 

5A 2 = -A- + ^-h {2) ) sin((fc - 1)6). (26) 



4. Symmetries and Vortex Scattering 

In the Abelian Higgs model, the zero modes can be used to study vortex 
scattering in the context of the slow-motion approximation ||10|| . The idea of 
this approximation is that for low velocities, at each point in time the fields, 
Ai, A 2 , and <f), are given by one of the static solutions; i.e., as an ansatz for 
these functions one can choose the family of static solutions after making 
the parameters time dependent. A is then determined from its equation of 
motion, which is of zero order in t. Finally, with this ansatz the action is 
minimized. 

For our purpose, this idea can be implemented as follows. Near the 
rotationally symmetric vortices, we expand the fields, 

f = ft + s\t)5(t)\ A' = A l + s 2+l (t)5^ (27) 

Here, <p = 4> l + ift 2 and A 1 are the the static solutions (7) and (8). 5(f) 1 and 
5A 1 are the fluctuations (21,22,23) or (24,25,26). The equation, which has to 
be solved for Aq, is the first (for v = 0) of the following three equations: 

dM -\<P\ 2 Y P ~\(V 2 -\<P?)F» U + i(p-l)D^<f>*D^<f>)) + (p-l) 

x( ^_| |2 )2p -4 (((r?2 _ |0|2)F „ + ,( p _l) £) lM^]0)(^ £) ^ + ^*)) 

-ip(2p - l)(ri 2 - \4>\ 2 ) 2p - 2 x {4>D V 4>* - 4>*D v (t>) = 0. (28) 

(The other two equations (for v = 1,2) are the second-order equations of 
motion for A\ and A 2 , which we will need later.) 
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The dynamics, as discribed by the functions s a (t), is now given by the 
Lagrange function, 

lw = / 101 2 ) 2(p - 2) (2((^ 2 -i0i 2 )^ +.(p-iw A]0)% onN 

+4p(2p-l)(r / 2 -|0| 2 ) 2 | J D o 0| 2 ) d 2 x. 

The functions s a are found by solving the Euler-Lagrange equations of this 
Lagrange function. Obviously, the slow-motion approximation still leaves us 
with very complicated equations to solve. With our ansatz (27), however, we 
only attempt to study the neighbourhood of rotationally symmetric vortices; 
i.e., we can neglegt higher order terms in the functions s a . Even though we 
are allowed to linearize in s a , we were not able to solve Eq. (28) for v — 0. 

In the Abelian Higgs model, the Gauss equation corresponds to Eq. (28) 
for v = 0. Following the same steps we have just discussed for our models, in 
the Abelian Higgs model one finds that A = 0, and that the s a are functions 
linear in time, in the neighbourhood of the rotationally symmetric vortices. 
For the modes with real smooth 5(f) this implies - scattering, if only times 
shortly before and shortly after the vortices coincide are considered. We 
see that the slow-motion approximation yields more results in the Abelian 
Higgs model. Another difference between the models discussed here and the 
Abelian Higgs model is that in the Abelian Higgs model the validity of the 
slow-motion approximation has been proved rigorously ||TT |. 



Instead of pursuing the slow-motion approximation any further, we will 
now formulate a corresponding Cauchy problem and use the symmetries to 
study the full Euler-Lagrange equations of the Lagrangian (1). We work 
with the vector, ip T = (A , d t A , Ai, d t A ly A 2 , d t A 2 , 0, d t 4>), and impose the 
following initial conditions: 

^(0, xf = (i , 0, At, 6A U A 2 , SA 2 , 0, 50) (30) 

Here, and A^ are the static solutions (7) and (8). 5<p and 5Ai are the 
fluctuation (21,22,23) for k = n. We concentrate on this type of fluctuations 
for simplicity, and because they lead to the interesting ^ symmetry in the n 
vortex scattering process. Aq is the solution of Eq. (28) for v — at time 
t = 0. 

The equations of motion we consider are the following second-order equa- 
tions: The second-order equation for A Q is, d tt A°+did t A l = 0. (This equation 
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follows from the Lorentz condition d^A^ 1 = 0.) The second-order equations 
for A\ and A2 are given by Eq. (28) for v — 1,2. The second-order equation 
for is its Euler-Lagrange equation, 



zd,(( V 2 - |0| 2 ) 2p -W - |0| 2 )i^ + i(p-l)I?^*I^](p-l)A^* 

+p(2p-l)(r ? 2 -|0| 2 ) 2 W) + i(p-2)(r / 2 -|0| 2 ) 2 ^ 3 x 

([(^-i^H^ + ^-i)^*^^] 2 + Ap{2p-i){r, 2 -\<p\ 2 f\D,<p\ 2 

+2(2p-l) 2 (r / 2 -|0| 2 )V " \(v 2 -m 2ip - 2 \((v 2 -\<P\ 2 ) 

xF^+i(p-l)D^<f>*D^<f>) x (F^-ip-l^D^*) 

-2 P {2 V -l){rf-\cj>f)\D^ ( j ) * + zp(2p-l)(r ? 2 -|0| 2 ) 2 ^ J D^ 

-2(2p-l) 2 (7 7 2 -|0| 2 ) 3 0*) = 0. (31) 



The Lorentz condition and Eq. (28) for v = are considered as constraints, 
which, by our choice of initial data, are satisfied at t = 0. In the Abelian 
Higgs model, it has been proved that the analogous Cauchy problem has a 
unique solution, and that the constraints are propagated f7| . In the following, 
we will assume that also in the present case a unique solution exists, although 
a rigorous proof is still missing. (Because of the complexity of the equations, 
to rigorously prove existence of a unique solution seems a very difficult task.) 

For a given solution ip(t, x), we define the functions, ipi(t, x) = Miijj(t, X(i)) 
for i — 1, 2, where xm = Sx with 




(32) 




\ 
ABO 
-B A 
1/ 



{xi, -x 2 ) 



. The matrices Mj are defined as, 



/ 



( 



1 
1 



) 



A = cos — /, B = sin — 

n n 



/, (33) 



and 



M 2 



I -I \ 
0-/00 

0/0 

V 00c/ 



CV = V* . 



(34) 
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We can now show that, if ip(t,x) is a solution of the Cauchy problem, so 
are ipi(t, x) and ^(t, x). (To find the symmetries of the initial data A o (0, x), 
we have assumed that Eq. (28) for v = has a unique smooth solution 
with asymptotic decay sufficient to satisfy the finite-energy condition.) The 
uniqueness of the solution of the Cauchy problem now implies that, actually, 
■0, ifti, and ip2 are all the same function. This, in turn, implies that functions 
like \4>\ 2 , Ffj, or the energy density S are invariant under a ^ rotation, and 
under a reflection w.r.t. the Xi-axis. This leads to the following conclusion: 
If by using functions like |0| 2 ,F?- or S, there is a way of defining the positions 
(xi(t), x^{t)), a = l,...,n, of exactly n separate vortices, these n positions 
must lie on n radial lines separated by an angle — with equal distance from 
the origin. (As in Ref. 7, we can use the the minima of |$| 2 to define these 
positions, near the rotationally symmetric vortices.) Furthermore, one of 
these radial lines must be the positive Xi-axis, or make an angle - with the 
positive Xi-axis. Any vortex that does not satisfy these conditions immedi- 
ately leads to 2n — 1 other vortices, because of the symmetries of our solution. 
For continuous solutions, these positions will change continuously such that 
at t — the n positions coincide, and after the collision the vortices move 
again on the radial lines just described. Therefore, they can either go back 
on the radial lines they came in on, or go back on radial lines shifted by an 
angle -. We will study a further symmetry to show that the second case is 
realised. 

The last transformation we study is x — > Mx, where M is the orthogonal 
matrix 




(35) 



Under this transformation the initial data change as follows: 



ip(0,Mx) = M 3 V(0,f), 



with 




C 
D 





—D 

C 




\ 




M 3 






V o 






-a J 



, C = cos —a , D = sin —a. 

n n 



(36) 
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(We have again assumed the uniqueness of the smooth finite-energy solution, 
A)(0, x), of Eq. (28) for v = 0.) From the uniqueness of the solution of 
the Cauchy problem, ip(—t,Mx) = M 3 vp(t,x) follows, and we see that the 
functions, |0| 2 , Ff- and S, are invariant under the transformation (t,x) — > 
(—t, Mf). This establishes ^ scattering for n vortices. 

5. Conclusions 

The study of vortices in a hierarchy of generalized Abelian Higgs models, 
and the comparison with the Abelian Higgs model was continued. In previ- 
ous studies we had seen that the rotationally symmetric generalized vortices 
are qualitatively similar to, but quantitatively different from the vortices. 
Here similar structures were found in the neighbourhood of the rotationally 
symmetric vortices. We found that the zero modes of the generalized vor- 
tices have the same angular dependence as, but radial behaviour different 
from that of the zero modes of the vortices in the Abelian Higgs model. 

The slow- motion approximation turned out to be of limited value. (Since 
the time-independent solutions are not known explicitly even in the Abelian 
Higgs model, the slow-motion approximation is not very successful in this 
model either.) On the reasonable assumption that a certain Cauchy problem 
has a unique solution, we were, however, able to study the symmetries of 
certain solutions. Each solution describes a process where n vortices approach 
and form one structure, namely a rotationally symmetric vortex. From this 
structure n vortices emerge. The pattern the vortices create for time t is the 
same as that for time —t, after a - rotation. 
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